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Ongoing experimental eflorts to detect cosmic sources of high energy neutrinos are guided by the 
expectation that astrophysical accelerators of cosmic ray protons would also generate neutrinos 
through interactions with ambient matter and/or photons. However there will be a reduction in the 
predicted neutrino flux if cosmic ray sources accelerate not only protons but also significant num- 
bers of heavier nuclei, as is indicated by recent air shower data. We consider plausible extragalactic 
sources such as active galactic nuclei, gamma-ray bursts and starburst galaxies and demand consis- 
tency with the observed cosmic ray composition and energy spectrum at Earth after allowing for 
propagation through intergalactic radiation fields. This allows us to calculate the expected neutrino 
fluxes from the sources, normalized to the observed cosmic ray spectrum. We find that the likely 
signals are still within reach of next generation neutrino telescopes such as IceCube. 
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I. INTRODUCTION 

It has long been recognised that high energy pro- 
tons produced in cosmic ray accelerators would also gen- 
erate an observable flux of cosmic neutrinos, mainly 
through pion production in collisions with the ambi- 
ent gas or with radiation fields Neutrino telescopes 
such as AMANDA/IceCube 0, ANTARES Q, RICE [1], 
ANITA [5] and the Pierre Auger Observatory have 
presented initial results and are approaching the level of 
sensitivity thought to be required to detect the first high- 
energy cosmic neutrinos 01 ■ 

Candidate sources of high-energy neutrinos are of both 
Galactic and extragalactic varieties, the latter being ex- 
pected to dominate at the highest energies just as for 
the parent cosmic rays. Likely Galactic sources include 
microquasars Q and supernova remnants , while pos- 
sible extragalactic sources include active galactic nuclei 
(AGN) J^], gamma-ray bursts (GRBs) 11], and star- 
burst galaxies [l^ . If however the cosmic ray sources 
accelerate heavy nuclei, the resulting high energy neu- 
trino spectrum may be altered considerably from the all- 
proton picture which is usually assumed. Nuclei under- 
going acceleration can interact with radiation fields in 
or near the cosmic ray engine, causing them to photo- 
disintegrate into their constituent nucleons which can 
then proceed to generate neutrinos through photo-pion 
interactions. Hence if most of the accelerated nuclei are 
broken up into nucleons before they can escape from 
their sources, the neutrino spectrum will not differ much 
from that predicted for proton accelerators. However 
if the radiation fields surrounding the sources are not 
sufficiently dense to fully disintegrate cosmic ray nu- 
clei, fewer nucleons will be freed, leading to a reduced 
neutrino flux. In fact heavy nuclei can directly photo- 



produce pions on radiation fields but since the photo- 
pion production threshold is much higher than typical 
photo-disintegration thresholds, such interactions will be 
unimportant except at very high energies, well beyond 
the Galactic/extragalactic transition in the cosmic ray 
spectrum. Of course both nuclei and nucleons can scatter 
inelastically with ambient gas surrounding the sources to 
produce pions which subsequently decay into neutrinos. 

In this article, we explore the impact of primary nuclei 
on the generation of high energy neutrinos in plausible 
extragalactic sources of cosmic rays. In Sec. [Hi we re- 
visit the "Waxman-Bahcall bound" on the high-energy 
cosmic neutrino spectrum and generalize their argument 
to include cosmic ray nuclei. In Sec. IIIII we estimate 
the photo-nuclear energy loss lengths for three suggested 
high-energy cosmic ray accelerators: AGN, GRBs, and 
starburst galaxies. We find that from AGN, nuclei with 
energies below ^ 10^^ eV can escape largely intact, while 
by 10^^'^ eV most iron nuclei will suffer violently and 
only nucleons or light nuclei are expected to escape. In 
the case of GRB, most energetic nuclei will undergo com- 
plete disintegration and only nucleons can be emitted 
as cosmic rays. Starburst galaxies by contrast are ideal 
candidates for the emission of ultra-high energy cosmic 
ray nuclei, as very few nucleons are dissociated. (Re- 
cently, accretion shocks around galaxy clusters have also 
been invoked in this context In Sec. IIVI we dis- 

cuss the implications of these results and derive the neu- 
trino fluxes associated with the observed ultra-high en- 
ergy cosmic rays under the assumption that none of the 
sources are hidden. To normalize the fluxes we require 
that their energy spectrum and chemical composition at 
Earth, after propagation through intergalactic radiation 
backgrounds, are consistent with experimental data from 
Auger and HiRes. In Sec.|V]we show that the expected 
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cosmic neutrino fluxes are still detectable by kilometer 
scale neutrino telescopes, such as IceCube. Our conclu- 
sions are presented in Sec. |Vll 



II. THE WAXMAN-BAHCALL BOUND 

Waxman and Bahcall pointed out some time ago 
that an upper bound can be placed on the diffuse flux 
of cosmic neutrinos assuming that they are generated in 
cosmologically distributed, optically-thin proton acceler- 
ators (see also Refs. [l^ and [l^)- In this Section, we 
discuss how this bound is affected if the accelerated par- 
ticles are nuclei rather than protons. 

First we briefly outline the original argument as ap- 
plied to proton accelerators. If the observed flux of ultra- 
high energy cosmic rays is the result of cosmologically 
distributed sources, then the energy injection rate in the 
10^^ — 10^^ eV) energy range can be inferred to be [l7j : 



CR 



CR 



dEcK 



.[lol^lo"l] 
ln(102VlOi9) 



10'^'' ergMpc^Vr"^ 



(1) 

where an energy spectrum oc -B ^ has been assumed and 
Eq = 10^^ eV. The energy density of neutrinos produced 
through photo-pion interactions of these protons can be 
directly tied to the injection rate of cosmic rays: 
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where tu is the Hubble time and e^r is the fraction of 
the energy which is injected in protons lost to photo- 
pion interactions. (The factor of 3/8 comes from the 
fact that, close to threshold, roughly half the pious pro- 
duced are neutral, thus not generating neutrinos, and 
one quarter of the energy of charged pion decays - 
TT"*" — > ^i^v^i — > e^Vf,v^v^- going to electrons rather than 
neutrinos.) Thus the expected neutrino flux is: 
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where the parameter accounts for the effects of source 
evolution with redshift, and is expected to be of order 
unity. The "Waxman-Bahcall bound" is defined by the 
condition e^r = 1- For interactions with the ambient gas 
(i.e., pp rather than p^ collisions), the average fraction 
of the total pion energy carried by charged pions is 2/3, 
compared to 1/2 in the photo-pion channel. In this case, 
the upper bound given in Eq. ^ is enhanced by 33% [3 • 
At production, the neutrino flux consists of equal frac- 
tions of Ve, v^j, and i>^. Originally, the Waxman-Bahcall 
bound was presented for the sum of and (neglect- 
ing Uf.), motivated by the fact that only muon neutri- 
nos are detectable as track events in neutrino telescopes. 



Since oscillations in the neutrino sector mix the differ- 
ent species, we chose instead to discuss the sum of all 
neutrino flavors. When the effects of oscillations are ac- 
counted for, nearly equal numbers of the three neutrino 
flavors are expected at Earth. Note that IceCube will be 
capable of detecting and discriminating between all three 
flavors of neutrinos . 

Electron antineutrinos can also be produced through 
neutron /3-decay. This contribution, however, turns out 
to be negligible. To obtain an estimate, we sum over the 
neutron-emitting sources out to the edge of the universe 
at a distance ^ 1 / Hq : 
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where dNn/dEn is the neutron volume emissivity and mn 
the neutron mass. Here, we have assumed that the neu- 
trino is produced monoenergetically in the neutron rest 
frame, i.e., eq ^ Sm{l — ml/S^)/2 ^ 0.55 MeV, where 
Sm ~ 1.30 MeV is the neutron-proton mass difference. 
An upper limit can be placed on dNn/dEn by assuming 
an extreme situation in which all the cosmic rays escap- 
ing the source are neutrons, i.e.. 
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With the production rate of ultra-high energy protons 
^[10", 10=1] ^ 5 X iQiA gj.g Mpc-3 yr-i [l|], and an as- 
sumed injection spectrum dNn/dEn oc i?,7^, Eq. (jH) gives 
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which is about three orders of magnitude below the 
Waxman-Bahcall bound in Eq.Q. 

If the injected cosmic rays include nuclei heavier than 
protons, then the neutrino flux expected from the cosmic 
ray sources may be modified. Nuclei undergoing acceler- 
ation can produce pions, just as protons do, through in- 
teractions with the ambient gas, so the Waxman-Bahcall 
argument would be unchanged in this case. However if 
interactions with radiation fields dominate over interac- 
tions with matter, the neutrino flux would be suppressed 
if the cosmic rays are heavy nuclei. This is because the 
photo-disintegration of nuclei dominates over pion pro- 
duction at all but the very highest energies. Defining n 
as the fraction of nuclei heavier than protons in the ob- 
served cosmic ray spectrum, the resulting neutrino flux 
is then given by: 
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The Waxman-Bahcall bound can be turned into a flux 
prediction by making further assumptions. For example 
one can assume that all charged particles remain trapped 
within the acceleration region and only neutrons are able 
to escape. If this were the case, the energy fraction of 
cosmic ray protons lost to photo-pion interactions can 
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be easily obtained from single pp or p7 collisions, yield- 
ing e^r ^ 0.2 — 0.6. This flux estimate, however, does 
not account for high energy proton production at the 
source edge (which would have a large probability of es- 
caping the acceleration region), threshold effects (near 
the photo-pion production threshold, the assumed re- 
lationship between the average energy of the incoming 
proton and that of the outgoing neutrino can be signif- 
cantly altered), and energy degradation of the charged pi- 
ous propagating through the (possibly strong) magnetic 
fields in the plasma. In the remainder of this paper, we 
will estimate the flux of cosmic neutrinos from specific 
astrophysical sites, taking into account all of these con- 
siderations and assuming that high energy cosmic rays 
are constituted of both protons and nuclei, in conformity 
with observational data on their spectrum and composi- 
tion. 



III. PHOTO-REACTION RATES FOR 
PROTONS AND NUCLEI: THREE CASE 
EXAMPLES 

It is likely that the bulk of the cosmic radiation is cre- 
ated as a result of some general magneto-hydrodynamic 
process which channels kinetic or magnetic energy of 
cosmic plasmas into charged particles. The details of 
the acceleration process and the maximum attainable 
energy depend on the time scale over which particles 
are able to interact with the plasma. Sometimes the 
acceleration region itself only exists for a limited pe- 
riod of time (e.g. supernovae shock waves dissipate af- 
ter about IC* yr). If the plasma disturbances persist 
for long periods, the maximum energy may be limited 
by the likelihood of escape from the region. If one in- 
cludes the effect of the characteristic velocity, /3c, of the 
magnetic scattering centres, the above argument leads 
to the so-called "Hillas criterion" for the maximum en- 
ergy acquired by a charged particle moving in a medium 
with magnetic field B, E,^g.^ ~ 2(3 c Ze B r^^, where tl ~ 
Z-^{B/fiGy^ eV) k£C is the Larmor radius 

of cosmic rays with charge Ze [21j. In what follows, we 
consider the radiation fields associated with three classes 
of potential cosmic ray sources and determine the parti- 
cle's energy loss rate, assuming that within these sources 
the trapping condition for efficient acceleration is fulfilled 
up to the highest energies 0. 

There are basically four ways in which ultra-high en- 
ergy cosmic rays interact with ambient photons: Comp- 
ton interactions, photo-pair production in the field of 
the nucleon/nucleus, photo-disintegration of the nucleus, 
and photo-production of hadrons (mainly pions). Al- 
though Compton scattering has no energy threshold, it 
results in a negligibly small energy loss rate for high 
energy cosmic rays, so we will not consider it further. 
Photo-pair and photo-pion production occur at center- 



relative contributions of these two processes to the to- 
tal energy losses is dictated by the ratio of the product 
of the corresponding inelasticity and cross-section. For 
a relativistic nucleon, this product has an average value 
of ^ 5 X 10~'^^cm^ from threshold up to 500 MeV in 
the case of photo-pair production, and an average value 
of 7 X lO'^'^cm^ from threshold up to ~ 200 GeV for 
the case of photo-pion production [23|. For nuclei, the 
energy loss rate due to photo-pair production is Z^/A 
times higher than for a proton of the same Lorentz fac- 
tor (2JI , whereas the energy loss rate due to photo-meson 
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production remains roughly the same (because the cross- 
section for photo-meson production by nuclei is propor- 
tional to the mass number A [2^ , while the inelasticity is 
proportional to 1/A). However, it is photo-disintegration 
rather than photo-pair and photo-meson production that 
determines the energetics of ultra-high energy cosmic ray 
nuclei p6j . Experimental data on photo- nuclear inter- 
actions are consistent with a two-step process: photo- 
absorption by the nucleus to form a compound state, fol- 
lowed by a statistical decay process involving the emis- 
sion of one or more nucleons. Details of how the relevant 
cross sections for all these processes are calculated and 
the energy loss rates implemented in our simulation can 
be found in Ref. [H. 



A. Active Galactic Nuclei 

Among known non-thermal sources in the universe, 
radio-loud AGN seem to be the most energetic [28] hence 
they have long been suspected to be likely accelerators of 
ultra-high energy cosmic rays [2^. At radio frequencies 
where very large baseline interferometers can resolve the 
emission regions on milliarcsecond scales, many radio- 
loud AGN exhibit highly coUimated jets of relativistic 
plasma with opening angles of a few degrees or less. 
AGN come in various guises according to the orientation 
of their radio jets and the characteristics of the circum- 
nuclear matter in their host galaxies. The most extreme 
versions are Fanaroff Riley radio-galaxies whose radio jet 
axes are almost in the plane of the sky, and blazars which 
have the radio jet axes aligned close to the line-of-sight 
to the observer (yielding a significant flux enhancement 
through Doppler boosting). 

A total of 66 blazars have been detected to date as GeV 
7-ray sources by EGRET [sO]. In addition, observations 
from ground-based air Cerenkov telescopes indicate that 
in more than ten of these sources the 7-ray spectrum ex- 
tends above a TeV [3l|. The non-thermal emission of 
these powerful objects shows a triple-peak structure in 
the overall spectral energy distribution. The first com- 
ponent (from radio to X-rays) is generally interpreted 
as being due to synchrotron radiation from a population 
of non-thermal electrons, the second (UV) component as 
blackbody radiation originating from the accretion disk, 
and the third component (7-rays) is explained either as 
due to inverse Compton scattering of the same electron 



4 



LU 
T3 



LU 




-2 2 4 

loQio Ey [eV] 

FIG. 1: The photon energy spectrum of a flaring blazar in the 
observer's frame. The normahzatiori is obtained through the 
relation, U^/E'^"''^ = E'^'"'^dN/AE^i\^^^^ (N.B.- such a nor- 
malisation does not give the photon number density actually 
measured in the observer's frame). 



population on the various photon fields traversed by the 
jet [34I, or due to development of electromagnetic cas- 



cades triggered by secondary pions (produced when the 
highly relativistic baryonic outflow collides with the pho- 
ton background fssj or diffuse gas targets moving across 
the jet [34]) that cool instanstaneously via synchrotron 
radiation. The characteristic single photon energy in syn- 
chrotron radiation emitted by an electron of energy E is 
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E~'^-^, and a 20, 000 K blackbody radiation from the ac- 
cretion disk. It is assumed that 10% of the Eddington 



luminosity is radiated as blackbody radiation and 1% as 
synchrotron radiation (37| . 

The cosmic ray acceleration process is assumed to oc- 
cur within a relativistic blob of plasma moving along the 
jet with Lorentz factor F ^ 10^'^. In the rest frame of the 
plasma, the blob is assumed to be spherical with radius 
R' = cAt' = TcAt, where At w 10* s indicates the typi- 
cal duration of a flare. (More exactly we should use the 
Doppler factor and average over the angle between the 
particle velocity and the boost velocity but ignoring such 
subtleties does not change our estimates substantially.) 
The energy density injected into such a blob can be esti- 
mated from the apparent bolometric luminosity, L-y , and 
the duration of flare 
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The spectrum of background photons is shown in Fig. [T] 
Note that since the blob is relativistic, L'^At' = L^At/T, 
implying a significant dilution of the energy density in the 
plasma rest frame to that expected in a non-relativistic 
scenario, 
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The fraction of the total energy in cosmic ray protons, 
E'p, expected to be lost within the acceleration region to 
pion production is just. 



where, B^q is the magnetic field in units of /iG and 
Eig = E/10^^ eV. For a proton of the same energy, this 
number is {mp/nief' ^ 6x 10^ times smaller. It is evident where 
that high energy gamma ray production through proton 
synchrotron radiation requires very large 0(100 G) mag- 
netic fields [IHl, and therefore we do not consider this 
process in our calculations. 

For the background radiation spectrum in AGN, we 
adopt the Bottcher blazar flaring state model [1^. The 
two relevant components are synchotron radiation peak- 
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(12) 



is the proton attenuation length due to interactions with 
the radiation field, cta is the cross-section, and K{E'p) is 
the inelasticity of a single collision. (Recall that primed 
quantities refer to values measured in the plasma rest 
frame). For protons interacting via the A-resonance we 
find. 
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where the inelasticity, Kp{Ep) « 0.2, is kinematically de- termined by requiring equal boosts for the decay products 
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FIG. 2: Comparison of the energy loss lengths for photo- 
pion production by protons, photo-pion production by iron 
nuclei and photo-disintegration of iron nuclei, as a func- 
tion of the proton/iron nucleus energy in the rest frame of 
the AGN plasma. The source size is normalized such that 
E'dx'/dE' = Kp{E'p)/eTr at the peak energy, which corre- 
sponds to a proton energy E' = E'j, 10^^'^ eV. The photo- 
disintegration and photo-pion production loss lengths for iron 
nuclei are normalized using the number density of photons 
given in Fig. [1] reduced by 1 /F^ to obtain the value in the 
plasma rest frame. 



of the A [3g| . A straightforward calculation shows that 
the required proton energy for A-resonant production at 
the peak of the photon distribution is, 
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where mA — 1232 MeV is the mass of the resonance and 
•nip the proton mass. In Fig. [5] we show the expected 
photo-pion production and photo-disintegration energy 
loss lengths for protons and nuclei. It is clear that the 
radiation field becomes thick to photo-pion interactions 
at an energy (in the plasma frame) of E'^ k, 10^^ eV 
(which corresponds to Ep « 10^^'^ eV in the observer's 
frame) , whereas it becomes opaque to the propagation of 
nuclei about a decade lower in energy. One can also ver- 
ify that photo-disintegration dominates over photo-pion 
production by nuclei at all energies, hence neutrino pro- 
duction through photo-nuclear interactions can be safely 
neglected. 



B. Gamma Ray Bursts 

GRBs are flashes of high energy radiation that can 
be brighter, during their brief existence, than any other 
source in the sky. The bursts present an amazing va- 
riety of temporal profiles, spectra, and timescales [s^ . 
Our insights into this phenomenon have been increased 



dramatically by BATSE observations of over 2000 GRBs, 
and more recently, by data from SWIFT [i^l . 

There are several classes of bursts, from single-peaked 
events, including the fast rise and exponential decay- 
ing (FREDs) and their inverse (anti-FREDs), to chaotic 
structures [4l|] . There are well separated episodes of emis- 
sion, as well as bursts with extremely complex profiles. 
Most of the bursts are time asymmetric but some are 
symmetric. Burst timescales range from about 30 ms to 
several minutes. 

The GRB angular distribution appears to be isotropic, 
suggesting a cosmological origin [42]. Furthermore, the 
detection of "afterglows" — delayed low energy (radio to 
X-ray) emission — from GRBs has confirmed this via the 
redshift determination of several GRB host-galaxies [i^ . 

If the sources are so distant, the energy necessary to 
produce the observed events by an intrinsic mechanism 
is astonishing: about 10^^ erg of gamma rays must be re- 
leased in less than 1 second. The most popular interpre- 
tation of the GRB-phenomenology is that the observable 
effects are due to the dissipation of the kinetic en ergy 
of a relativistic expanding plasma wind, a "fireball" [44| . 
Although the primal cause of these events is not fully un- 
derstood, it is generally believed to be associated with the 
core collapse of massive stars (in the case of long duration 
GRBs) and stellar collapse induced through accretion or 
a merger (short duration GRBs) [45i]. 

The very short timescale observed in the light curves 
indicate an extreme compactness that implies a source 
which is initially opaque (because of 77 pair creation) 
to 7-rays. The radiation pressure on the optically thick 
source drives relativistic expansion, converting internal 
energy into the kinetic energy of the inflating shell. Bary- 
onic pollution in this expanding flow can trap the radi- 
ation until most of the initial energy has gone into bulk 
motion with Lorentz factors of F ~ 10^ — 10'^ [i^. (In our 
calculations we set F = lO'^ ''^). The kinetic energy can be 
partially converted into heat when the shell collides with 
the interstellar medium or when shocks within the ex- 
panding source collide with one another. The randomized 
energy can then be radiated by synchrotron radiation and 
inverse Compton scattering yielding non-thermal bursts 
with timescales of seconds. Charged particles may be 
efficiently accelerated to ultra-high energies in the fire- 
ball's internal shocks, hence GRBs are often considered 
as potential sources of cosmic rays [13] . 

To describe the radiation fields associated with GRBs, 
we adopt a standard broken power-law spectrum: 
AN^/dE^ oc E~^ , where (3 —\,2 respectively at energies 
below and above the break energy, iJ^rcak _ ^ MgV [ll[, 
which fits the BATSE data well 48]. In many ways the 
situation is similar to the blob of emitted plasma in the 
AGN model. However, in the fireball's comoving frame, 
a spherical shock expands relativistically in all directions 
(with Lorentz factor F), and thus a change in geome- 
try is required. In a frame moving with Lorentz fac- 
tor F towards the observer, the shock has thickness AR' 
(=i?'/F), where R' is the initial size of the compact ob- 
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ject before the fireball phase. In the observer's frame, the energy density, (7^, of a shell of luminosity, L^^, radius 
shock is further compressed into a thin shell of thickness R' , and thickness Ai?' is fi^]. 
AE'/r (=i?7r2). Therefore, the fiducial value for the 
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where i?' = cAi and Ai?' = F cAt. The spectrum of this radiation field is shown in Fig. O 
The fraction of energy deposited in the GRB by an ultra-relativistic proton of energy E'^ is, 

For protons interacting through the A-resonance with photons of energy E'^^'''^^ ~ lOkeV, Eq. yields: 
f 0.2 \ r L^At\ r CTA \ n 9 /^"^^^'^V '^'^ \ /I MeV\ /2 X 10-3 s\ 



In Fig. 13] we show the expected photo-pion production 
and photo-disintegration energy loss length for protons 
and nuclei. This figure highlights that (with our as- 
sumptions) GRBs are expected to be optically thin to 
P7 interactions below 10^^ eV, with protons undergo- 
ing at most one interaction before leaving the source re- 
gion. It is worthwhile pointing out that although the 
nucleus photo-disintegration and photo-pion production 
loss lengths become comparable at laboratory (observer) 
energies E > 10^^ eV, most of the charged pions pro- 
duced in this energy regime readily lose energy through 
synchroton radiation in the strong i3-field of the fireball. 
For example, at 10^^ eV only 1% of the produced pions 
decay before losing a significant fraction of their energy. 
Since the neutrino flux falls by more than two orders of 
magnitude per decade increase in energy, the contribu- 
tion from pions suffering energy degradation would be 
negligible, and therefore can be ignored. 



C. Star bursts 




-6 -4 -2 2 4 6 



loQio [eV] 

FIG. 3: The photon energy spectrum of a GRB in the ob- 
server's frame. The normalization is obtained from the con- 
dition U^/E^'"'''' ^ E^''"'''dN/dE^\^^^^^ (N.B.- such a nor- 
malisation does not give the photon number density actually 
measured in the observer's frame) 



Starbursts, galaxies undergoing an episode of large- 
scale star formation, have also been proposed as sources 
of ultra- high energy cosmic rays '50]. These environ- 
ments feature strong infrared emission by dust associated 
with high levels of interstellar extinction, strong UV spec- 
tra from the Lyman a emission of hot OB stars, and con- 
siderable radio emission produced by recent supernova 
remnants. The central regions of starburst galaxies can 
be orders of magnitude brighter than those of normal spi- 
ral galaxies. From such an active region, a galactic-scale 
superwind (driven by the collective effect of supernovae 
and winds from massive stars) can conceivably accelerate 



cosmic ray nuclei to ultra-high energies. 

The Lyman-a background is powered by the rich OB 
and red supergiant stellar populations within the inner 
core of these galaxies, which typically has a radius R ~ 
100 pc. The average density in the region depends on the 
temperatures Tqb and Tsg of O, B, and red supergiant 
stars, respectively, and the dilution due to the inverse 
square law. Specifically, for a region with Nqb OB stars 
and iVsG red supergiant stars, the photon density was 
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where 



e/T 



1 



(18) 



(19) 



is a Bose-Einstein distribution with energy e and temper- 
ature T (normalized so that the total number of photons 
in a box is J n{e)de). Here, i?oB(SG) is the radius of the 
OB (SG) stars and the factor 9/4 comes from averaging 
the inverse squares of the distance of an observer to uni- 
formly distributed sources in a spherical region of radius 
R We assume the OB and red supergiant stars to 
constitute 90% and 10%, respectively, of the total popu- 
lation of 27,000 [H. 

We also adopt a single component dust model in which 
clumpy dust surrounds the starburst region, in ther- 
mal equilibrium with it and heated to a temperature of 
30 K [54] . About 90% of the stellar light is assumed to be 
reprocessed into IR radiation and the resulting thermal 
photon spectrum is shown in Fig. [5] 

Iron nuclei with Lorentz factors up to ~ 10^ may be ac- 
celerated in supernova explosions 55]. Despite the star- 
burst region being only about ~ 100 pc in size, most of 
these nuclei may be trapped there through diffusion in 
milli-Gauss magnetic fields [H^ for ~ 10^ yr [53]. How- 
ever since magnetic rigidity increases with Lorentz factor, 
the diffusion time would decrease with rising energy. In 
our analysis we adopt the cosmic ray time delay 
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FIG. 6: Comparison of the photo-pion and photo- 
disintegration energy loss lengths as a function of the pro- 
ton/iron nucleus energy in a starburst galaxy. 



as expected due to trapping in a Kolmogorov spectrum 
of magnetic field fluctuations [ssj . 

Combining the above estimates of the propagation 
time, the size of the starburst region and the number 
density of UV and IR photons, we find the proton and 
nucleus photo-production energy loss lengths shown in 
Fig. [6l It is seen that the medium is (almost) trans- 
parent to the propagation of cosmic rays. Photo-pion 
production by nuclei will not be important until very 
high energies (Lorentz factors ^ 10^"), hence can be ne- 
glected. 

Before proceeding further, it is important to note that 
pions produced through hadronic interaction with the 
ambient gas would also contribute to the flux of neu- 



8 




oEpe [eV] 



emission of ultra-high energy cosmic ray nuclei, as very 
few nucleons are dissociated even at the highest energies. 

We can now calculate the neutrino flux associated with 
the observed ultra-high energy cosmic rays, taking into 
account observational data on their energy spectrum and 
showering characteristics which suggest that they consist 
of a mixture of protons and heavy nuclei^ rather than 
being purely protons. In other respects we make the same 
assumptions as Waxman and Bahcall [l^], e.g. that the 
extragalactic sources are not hidden. 



IV. THE COSMIC NEUTRINO FLUX FROM 
EXTRAGALCTIC SOURCES 
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FIG. 7; The top panel shows the average number of nucleons 
photo-dissociated from an iron nucleus of a given energy for 
the radiation field spectra of AGN, GRBs, and starbursts. 
The lower panel shows the starburst case in more detail. 



trinos. The exact intensity of this "pp hadronic compo- 
nent" is currently under debate 0, HH but an upper 
bound can certainly be placed from the non-observation 
of 7-rays from nearby starburst galaxies [60j . In this re- 
gard, the improved sensitivity of GLAST over EGRET is 
significant — detection of even a single starburst galaxy 
by GLAST would imply that such objects make a con- 
siderable contribution to the diffuse flux of extra-galactic 
neutrinos 61]. To be conservative in our calculations we 
do not consider contributions from the pp channel. 

Our essential results so far can he summarized as 
follows with reference to Fig. 

(1) In AGN, nuclei with energies below about 10^^ eV can 
escape largely intact, while more energetic nuclei undergo 
an increasing degree of disintegration; by 10^^'^ eV, most 
iron nuclei from AGN will suffer violently, and only 
nucleons or very light nuclei are expected to escape. 

(2) In GRB, most nuclei will undergo complete disinte- 
gration and only nucleons can be emitted as cosmic rays. 

(3) Starburst galaxies are ideal candidates for the 



As we have seen above, cosmic ray nuclei which escape 
from the acceleration region unscathed do not contribute 
to the source neutrino spectrum. Therefore, in order to 
estimate the latter, we first need to determine the frac- 
tion of heavy nuclei in the high energy cosmic rays arriv- 
ing at Earth, which we denote by n. 

This issue is closely connected with that of the "cross- 
over" energy at which a transition occurs between Galac- 
tic and extragalactic cosmic rays. It would be natural 
to expect a flattening of the energy spectrum at this 
point as the harder subdominant extragalactic compo- 
nent takes over from the softer Galactic component. Such 
an "ankle" was indeed observed by Fly's Eye at 
~ XQ18.5-18.7 gY .^^j^igj^ jg roughly the energy at which 

the (proton) Larmor radius begins to exceed the thick- 
ness of the Milky Way disk and one expects the Galactic 
component of the spectrum to die out. The end-point of 
the Galactic flux ought to be dominated by heavy nuclei, 
as these have a smaller Larmor radius for a given energy, 
and the data is indeed consistent with a transition from 
heavy nuclei to a lighter composition at the ankle. 

However recent HiRes data indicate that this change 
in the cosmic ray composition occurs at a much lower en- 
ergy of ^ 10^^'^ eV [6J|, where the spectral slope steep- 
ens from E"'^ to E~'^-^ [63|. This "second knee" in the 
spectrum, recognised originally in AGAS A data [6^1 , can 
be explained [66| as arising from energy losses of extra- 
galactic protons over cosmic distances due to e+e" pair- 
production on the cosmic microwave background (CMB). 
The "ankle" is now interpreted as the minimum in the 
e'^e~ energy-loss feature ^67i] and this requires that there 
be an almost total absence of iron nuclei 0.05% for 
E > 10^^ eV) in extragalactic cosmic rays w]. 

However the data on the composition inferred from 
different analyses of the characteristics of cosmic ray 
air showers do not support such a simple picture [69| . 
For example the depth of shower maximum Xmax at 
E > 10^^ eV is not consistent with a pure proton com- 
position, even allowing for the uncertainties in mod- 
elling hadronic interactions at such high enegies. This 
is seen in Fig. [5] where we show the evolution of Xmax 
for both proton and iron showers as obtained from ex- 
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FIG. 8: The top (bottom) panel shows the depth of shower 
maximum, Xmax, for cosmic rays with energy spectrum oc 
E-" which are a mixture of 90% (75%) protons and 10% 
(25%) iron. The solid, dashed, and dotted lines are the pre- 
dicted values of Jfmax for a pure iron reaching the Earth's 
atmosphere, pure proton, and the mixed composition (at the 
source), respectively. The spread in the predictions for dif- 
ferent hadronic interaction models is also shown. The ex- 
perimental data are from the Fly's Eye ^62], HiRes [t^I and 
Yakutsk [7^1 experiments. 



tensive air shower simulations using the program COR- 
SIKA (version 6.20) [t^ along with experimental data 
from Fly's Eye, HiRes and Yakutsk. The predictions 
of three different hadronic interaction models (DPM- 
JET [761], SIBYLL '{ITi], and QGSJET II ^) are shown 
so that an impression of the modelling uncertainties may 
be gained from the spread between them. 

In fact, recent analyses [13, [T^I show that the Xmax 
data, as well as the data on the energy spectrum from 
both HiRes [7l!| and the Pierre Auger Observatory [t^ . 
can simultaneously be reproduced if the extra-galactic 
cosmic rays contain a substantial fraction of heavy nuclei, 
and the Galactic/extra-galactic transition occurs, as was 
believed originally, at the "ankle" in the spectrum. 

In previous work (27j we undertook a detailed calcula- 
tion of the intergalactic propagation of ultra-high energy 



heavy nuclei through known cosmic radiation fields, in 
order to find the energy spectrum and composition at 
Earth. To keep things simple we now consider only pro- 
tons and iron nuclei. Armed with our previous results, 
we estimate the fraction of iron nuclei, k, in extragalac- 
tic cosmic rays, by finding the linear combination of pure 
iron and pure proton energy spectra that matches the 
Auger data [23 best, after propagation effects are ac- 
counted for. Our results, shown in Fig. |9l indicate that 
a mixed composition with 75% protons and 25% iron 
nuclei best reproduces the data (although the fit qual- 
ity is poor, suggesting correlated and/or overestimated 
uncertainties). As seen in Fig. |8l this particular mixed 
composition at source also allows a good match to the 
ATinax data. A better fit can perhaps be obtained by con- 
sidering a more complex composition (e.g. reflecting the 
composition of matter in plausible sources) but the qual- 
ity of present data is not good enough to warrant this. 
Even in this simple 2-component model, we obtain a dis- 
tinct improvement in the fit by increasing k from 10% to 
25%, and this is consistent with the energy spectrum. 

By performing a Monte Carlo simulation for cosmic ray 
protons and nuclei propagating through the sources con- 
sidered, the neutrino flux produced in each of the three 
case examples can now be calculated. The results shown 
in Fig. [TO] are obtained assuming a cosmological distri- 
bution of sources which accelerate cosmic rays with a 
spectrum oc Using the emissivity of AGN, the 

corresponding peak neutrino flux is, 

<P^iE^) = 2.3x10"^ (l-e-"'/^"^ 

X £;-2.4Qgy-lj,j^-2g-lg^-l ^ (^21) 

where Ei, is in GeV and we have normalized to a cosmic 
ray production rate of: 

4'r ' = 2.2 X 10^4 erg Mpc" ^r"' , (22) 

as indicated by Auger data. In the case of GRBs, because 
of the strong magnetic fields in the plasma, the neutri- 
nos are created by parent protons with energies below 
IQ-'^* eV. Following Ref. [69], we assume in this case a 
break in the proton injection spectrum (oc E^^ below 
10^* eV), so that the corresponding peak neutrino flux 
saturates the Waxman-Bahcall bound, yielding 

$^(£:^) = 1.0x10^^ (l-e-"'^^"^ 

X £:-2GeV"^cm"2g-igj.-i (23) 

Note that this normalization differs from that in Eq. ([3]) 
because of the different energy range of the source 
injection spectrum required to accomodate the Auger 
data [zil. 

We do not show the neutrino flux expected from star- 
burst galaxies because it is negligible compared to that 
from AGN or GRBs. The interaction length of protons in 
starburst galaxies is ~ 500 times the source size or larger 



10 



100 



10 



C/3 

'e 

o 

> 

o 

o 

LU 



0.1 



0.01 



100 




18.5 19 19.5 20 20.5 21 

loQio E [eV] 



75% p + 25% Fe 
P 
Fe 
Auger 
AE 




a = 2.4 



0.01 



18 18.5 19 19.5 20 20.5 21 

loQio E [eV] 

FIG. 9: Best fit of tfie Auger energy spectrum '72\, assuming a 
constant comoving density of sources which emit both protons 
and iron nuclei with energy spectrum oc E~" . In the fitting 
procedure we used data above iJ_ = 10^*'^ eV (top panel) 
and above E- = lO^* '^ eV (lower panel), finding x^/d.o.f. = 
52.5/15 and x^/d-o.f. = 41.6/13, respectively. The favoured 
mixture is 10% (25%) iron and 90% (75%) for E- = 10^** '^ eV 
{E- = 10^*'^ eV) — the dependence on E- arises because 
lower energy points have smaller statistical uncertainties. The 
horizontal error bars indicate the systematic uncertainties in 
energy measurement. 



at all energies (see Fig. [B]), leading to e^r < lO^'^, while 
AGN and GRBs have e,r ^ 1- Thus the uncertainties in 
the modeUing of starburst galaxies are immmaterial in 
this context. 

The required bi-modal composition at source (75% 
protons plus 25% of heavy nuclei) does however require 
that sources such as starburst galaxies accelerate most 
of the iron nuclei, with most of the protons coming from 
AGN and GRBs. In Fig. [Til wc show the sum of the 
neutrino fluxes produced by AGN and GRBs assuming 
they contribute approximately equally. Not surprisingly, 
the predicted diffuse flux in this simple model. 
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is just the expectation given in Eq.(I7]), which was based 
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FIG. 10: The neutrino spectra produced by protons and iron 
nuclei being accelerated in AGN (top) and GRBs (bottom). 
The fluxes have been normalized following the Waxman- 
Bahcall prescription, namely assuming that all ultra-high en- 
ergy particles observed on Earth are protons. 



on the assumption of complete trapping of charged par- 
ticles (recall that n = 0.25 and, for single P7 collisions, 
~ 0.2). Of course a more sophisticated estimate 
can be made using our results in Fig llOl when we know 
more about the relative contributions from the different 
possible sources to the overall cosmic ray flux, as well 
as the relative (possibly energy dependent) weighting of 
heavy nuclei with respect to protons. Being conserva- 
tive, presently we can only argue that the overall cosmic 
neutrino flux should be reduced by about 75%. 

In closing, we stress that the diffuse neutrino flux has 
an additional component originating in the energy losses 
of ultra-high energy cosmic rays en route to Earth. The 
main energy loss process here is photo-pion production 
in the CMB, which causes the steepening of the cosmic 
ray spectrum beyond 10^^-'' eV. The decay of charged 
pions produced in this process results in a diffuse flux of 
"cosmogenic" neutrinos 'so'], which is comparable to the 
fluxes shown in Fig. 1101 If there are heavy nuclei in ultra- 
high energy cosmic rays then they will preferentially lose 
energy through photo-disintegration rather than photo- 



11 



1e-06 




1e-12 



loQio E, [GeV] 



FIG. 11: The neutrino spectrum produced by protons un- 
dergoing acceleration in AGN and GRBs which are assumed 
to contribute equally to the cosmic ray protons observed at 
Earth. The Waxman-Bahcall bound, as obtained from Eq. Q 
by setting = = 1, is also shown for comparison. 



pion production, so the cosmogenic neutrino flux will also 
be suppressed as has been discussed elsewhere [HI, [s^ . 

Equipped with the fluxes shown in Figs. [TUl and [TTl we 
now proceed to determine the expected event rates in a 
next generation neutrino telscope such as IceCube. We 
do not consider the cosmogenic flux so these rates are 
conservative. 



V. EXPECTED EVENT RATES AT ICECUBE 

The IceCube neutrino telescope is currently under con- 
struction at the South Pole. When complete, it will com- 
prise a cubic-kilometer of ultra-clear ice about 2 km be- 
low the surface, instrumented with long strings of sensi- 
tive photon detectors [s^. We focus on this experiment 
as it is the most advanced of its kind, however the event 
rates estimated below will also apply to e.g. the KM3 un- 
dersea experiment being planned for the Mediterranean. 

IceCube is designed to observe muon tracks and show- 
ers produced by neutrino charged current (CC) and neu- 
tral current (NC) interactions in and around the instru- 
mented volume. The probability of detecting a neutrino 
passing through the detector from its muon track is given 
by 



Pu — ^zenith) 



zenith ), (25) 



where cr^MiEy) is the charged current neutrino- nucleon 
cross section [84|, n is the number density of nucleons in 
the ice, and the muon range R^{E^, 9 zenith) is the aver- 
age distance traveled by a muon of energy E^ before it is 
degraded below some threshold energy (taken to be 100 
GeV). This quantity depends on the zenith angle of the 
incoming neutrino as only quasi-horizontal or upgoing 
events can benefit from longer muon ranges. At the en- 



ergies we are most concerned with, the majority of muon 
events will be quasi-horizontal. 
The expected muon event rate is 



dt 



= j dE^ dn<l>^^{E^) P^^^,{E^,e^cnith) A^s ; 



(26) 



where {E,^) — ^,y{E^)/3 is the flux of muon neutrinos 
and Acs ~ 1 km^ is the effective area of the detector [sl] . 
Similarly, the expected number of shower events is 



rli 



dE„ dn $, 



.{E.)<j':T''\e,)VoS, (27) 



where a 



CC(NC) 
i/AT 



is the CC (NC) neutrino-nucleon cross 
section and Ves is the effective volume for detection of 
showers [ssj- The sum is over and CC interac- 
tions and all NC interactions, and we assume ^^^{E^) = 
^i^{E,y)/3 where a — e,^,T. Electron neutrino, CC in- 
duced showers carry all of the incoming neutrino energy, 
whereas muon neutrino (CC or NC) induced showers 
carry away an energy of (1 — y) S^, where y is the in- 
elasticity of the neutrino interaction. Tau neutrino NC 
showers have an energy (1 — y) E^, whereas tau neutri- 
nos CC events with energies below E^ ~ PeV generate 
showers with the full energy of the incoming neutrino. 
We assume that only showers with energies greater than 
3 TeV can be identified at IceCube. 

In Table 1, we show the predicted event rates for the 
various cosmic ray accelerators considered earlier. We 
also show the effect on the rates if the threshold for 
both muons and showers is taken to be 100 TeV — a 
cut at this higher energy is adequate to eliminate es- 
sentially all background from muon bremsstrahlung ra- 
diation near the detector and from muons produced in 
cosmic ray showers in the atmosphere. Moreover, the 
steeply falling flux of atmospheric neutrinos is negligible 
above this energy, so this cut generates a very pure sam- 
ple of extraterrestrial neutrinos. The labels "protons" 
and "iron" denote the species of particle which are accel- 
erated by the source (rather than what actually escapes). 

As expected, the event rates from GRBs are similar re- 
gardless of whether protons or iron nuclei are accelerated, 
as the latter are almost entirely disintegrated by the sur- 
rounding radiation fields. Values for starburst galaxies 
are not given in the table as the rates (from p'y interac- 
tions) are well below the sensitivity of IceCube and other 
next generation neutrino telescopes. 



VI. CONCLUSIONS 

We have studied the role that heavy nuclei play in the 
generation of neutrinos in possible astrophysical sources 
of high energy cosmic rays. During acceleration the 
nuclei may be completely photo-disintegrated into their 
constituent nucleons and we find this indeed happens in 
GRBs, resulting in the outgoing cosmic rays being pro- 
ton dominated. The neutrino flux is then left largely 
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TABLE L Predicted event rates at IceCube for various sources 
of high energy neutrinos, with the muon energy threshold 
set to 0.1 (100) TeV and the threshold for showers set to 
be 3 (100) TeV. The labels "protons" and "iron" denote the 
type of particle assumed to be accelerated, rather than what 
actually escapes (iron nuclei do not escape intact from GRBs, 
in particular) . The label "UHECR Best Fit" denotes the case 
where 75% of the ultra high energy cosmic rays arriving at 
Earth are protons (assumed to be accelerated equally by GRB 
and AGN) and 25% are iron nuclei (assumed to be accelerated 
by starburst galaxies). 



Source 


dMf,ldt 


dMs/dt 


AGN (protons) 


1.2 (0.34) yr^^ 


0.45 (0.089) yr"^ 


AGN (iron) 


0.23 (0.13) yr"^ 


0.045 (0.037) yr-i 


GRB (protons) 


16. (3.4) yr'i 


6.3 (2.2) yr-i 


GRB (iron) 


12. (2.8) yr-i 


4.5 (1.9) yr-i 


UHECR Best Fit 


6.5 (1.4) yr"^ 


2.5 (0.86) yr"^ 



The likely possibility of a substantial fraction of nu- 
clei in ultra-high energy cosmic rays implies therefore a 
somewhat reduced expectation for the neutrino flux from 
their cosmic sources. In particular, as the spectrum and 
elongation length in the atmosphere of ultra-high energy 
cosmic rays appears to be best fitted by a mixture of 
25% iron nuclei and 75% protons, the overall neutrino 
flux is reduced somewhat relative to the expectation for 
an all-proton cosmic ray spectrum. 

As next generation neutrino telescopes such as IceCube 
begin operation, it becomes increasingly important to re- 
fine the expectations for detection of high energy neu- 
trinos from the sources of cosmic rays. We have taken 
into account recent data on ultra-high energy cosmic rays 
to provide updated estimates of the neutrino fluxes to 
be expected from various types of possible extragalactic 
sources. The actual detection of cosmic neutrinos will in 
turn thus provide crucial information on the nature of 
the long sought sources of cosmic rays. 



unchanged from previous estimates which had ignored 
the possibility of nuclei being accelerated as well as pro- 
tons. At the other extreme, sources such as starburst 
galaxies hardly disintegrate accelerated nuclei, enabling 
such particles to escape and to contribute to the observed 
ultra-high energy cosmic ray spectrum, largely without 
contributing to the cosmic neutrino flux. In AGN the 
situation is in between, with nuclei being fully disinte- 
grated only at the highest energies, so the neutrino flux 
is suppressed at lower energies. 
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